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We analyze new spin effects in current-carrying state of superconductor-2D electron gas-su- 
perconductor (S-2DEG-S) device with spin-polarized nuclei in 2DEG region. The hyperfine in- 
teraction of 2D electrons with nuclear spins, described by the effective magnetic field B, produces 
Zeeman splitting of Andreev levels without orbital effects, that leads to the interference pattern of 
supercurrent oscillations over B. The spin-orbit effects in 2DEG cause strongly anisotropic depen- 
dence of the Josephson current on the direction of B, which may be used as a probe for the spin-orbit 
interaction intensity. Under certain conditions, the system reveals the properties of 7r-junction. 
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The spin-orbit (SO) and hyperfine (HF) interactions in 
GaAs heterojunctions and similar 2D quantum Hall sys- 
tems attract permanent theoretical and experimental at- 
tention. The hyperfine field of the nuclear spin subsystem 
acting upon the spins of charge carriers may reach 10 
G.EJ At low temperatures, tho-puclear spin relaxation time 
can be macroscopically longp so the nonequilibrium spin 
population in heterojunctioos, once created, is conserved 
during hundreds of seconds .El The Zeeman splitting com- 
bined with a strong spin-orbit coupling in GaAs/AlGaAs 
2DEG gives rise to a novel class of coherent phenomena, 
e.g., the spontaneous Aharonov-Bohm effecto 




FIG. 1. A model of the superconductor - 2DEG - supercon- 
ductor device based on GaAs/AlGaAs heterojunctionfl The 
normal n is directed towards the Al-doped layer. 

In this paper we discuss mesoscopic spin-orbit effects 
in Josephson current flowing across the S-2DEG-S struc- 
ture (Fig. 1) with polarized nuclei in 2DEG region. The 
transfer of the Josephson current through the normal 
conducting layer is provided by the Andreev reflection 
of electrons and holes at the NS-interfaces, which con- 
vert normal electron excitations into Cooper pairs in the 
superconducting banks. In a pure system with length d 



smaller than the electron scattering length, the interfer- 
ence between coherent electron states and retro-reflected 
hole states produces the set of spin-degenerate Andreev 
energy levels E\((fr), which depend on the quantum num- 
bers A and on the difference $ of thej-jOrder parameter 
phases in superconducting electrodes.u In short struc- 
tures, d <C £o (£o is the coherence length in the supercon- 
ductor), the Josephson current can be presented as the 
sum of currents transferred by individual Andreev bound 
states (see, e.g., Ref. Q and references therein), 

and must be sensitive to the HF and SO interaction which 
eliminate spin degeneracy of the Andreev levels. 

The contact HF interaction in a semiconductor is de- 
scribed by the HamiltonianH 



H h f = (8Tr/3)nBlh V IitrS(r - R t ). 

' * 1. 



(2) 



Here fxg is the Bohr magneton, 7^ is the nuclear magne- 
ton, and Ii, er, Ri, r are nuclear and charge carrier spins 
and positions, respectively. It follows from Eq. (2) that if 
the nuclear spins are polarized, (Y], Ii) 7^ 0, the charge 
carrier spins feel the effective HF field B which may cause 
spin splitting in 2J.1EG of the order of one tenth of the 
Fermi energy EpSsi The influence of the Zeeman split- 
ting solely on a supercurrent was studied first in Ref. O 
for the SFS junction (F denotes ferromagnetic metal). 
It was shown that the spin splitting suppresses the criti- 
cal current and produces its oscillations over the intrinsic 
magnetic field localized within the F-layer. 

The SO interaction of a charge carrier with the inter- 
face potential in GaAs/AlGaAs hetecew unctions is mod- 
elled by the Bychkov-Rashba termji30 



H so = (a/h) [crp] n, 



(3) 



where a = 0.6 x 10 9 eV cmj-Ap, hoies,ML3 and a = 
0.25 x 10 -9 eV cm for electrons ,t2lEJ er^, pi are the charge 
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carrier spin and momentum, n is the normal to the in- 
terface directed towards the Al-doped layer. 

The configuration proposed in Fig. 1 has the following 
characteristic features: 

i) the nuclear field B is localized outside the supercon- 
ductors and does not influence the pairing mechanism; 

ii) it affects only the electron spins and does not modify 
the space motion of charge carriers, whereas usual mag- 
netic field causes strong orbital effects and transforms 
Andreev levels into Landau bands ;EJ 

iii) since the SO interaction in Eq. (3) is spatially uni- 
form along the 2DEG plane, it does not suppress the 
supercurrent, in contrast to the SO scattering at nop-, 
magnetic impurities which acts as a depairing factor.E-51 
The situation changes in presence of the nuclear field. 

To calculate the Josephson current in Eq. (1), it is 
enough to find discrete eigenvalues of the BCS Hamilto- 
nian supplied by the HF and SO interaction terms, 



n = dv 



V' + (-B(p) + crA(r)^ + A(r)^+xjj+ + h.c. (4) 



where u is the vector of Pauli matrices, the spinor ip is 
composed from the annihilation operators ip* i °f ^ ne 
electron with a given spin, and E(p) = p 2 /2m— Ep is the 
energy of a normal electron excitation. For simplicity, we 
accept a step- wise model for the order parameter A(r), 



A(x) 



Aexp (i<j>signa;/2) , \x\ > d/2, 
0, |z| < d/2. 



(5) 



The vector A(r) in Eq. (4) describes the HF and SO 
interactions in accordance with Eqs. (2) and (3), 

k W = \h + w, \x\<d/2, h = ^B, w= |[pn]. (6) 

In order to consider uniformly the spin and electron- 
hole states of quasiparticles, it is convenient to express 
the Hamiltonian in Eq. (4) in terms of the 4-spinor field, 



Wf J 



</> + = (^t"> V'x") V'l) V't) 



H = - \ dV V +{r)H V {r), 



H = r z E{p) 



<j z [A(x)t + 
wcr + h z a z 



-A*(x)r-] 
+ T z h\\cr, 



\x\ > d/2, 
\x\ < d/2. 



(7) 



(8) 



(9) 



Here Tin is the component of Zeeman field in 2DEG plane; 
the Pauli matrices a act upon the spin states, while r- 
matrices operate in electron- hole space, e.g., 



G Z 

-a z 



(10) 



The problem of single-particle spectrum of the Hamil- 
tonian TL in Eq. (8) is equivalent to the solution of 
the Bogolyubov-de Gennes (BdG) equation for the 4- 
component wave function ^\(r) 



H^ x (r) = E x ^ x (r). 



(11) 



Assuming all characteristic energies A, h, w to be much 
smaller than Ep, we use the quasiclassical representation 
of *&\(r) as the product of rapidly oscillating exponent 
over the slowly varying envelope u(x): 

%>\(r) = exp(ispx + ip y y)u s (x), p = Jp\ - p 2 (12) 

where s = ±1 indicates two signs of x-component of the 
electron momentum. The spinor function u{x) obeys Eq. 
(11) with a reduced Hamiltonian, in which the quasi- 
classical approximation for the kinetic energy and SO 
operators is used, 

E{p) « svp x , w ps -r\pn], p = (sp,p y ,0), v = — . (13) 

The matching of the solutions of Eq. (11) at the NS 
interfaces, which are assumed to be completely transpar- 
ent, yields a dispersion relation, 



Ed 



= mi + arccos — 

hv A 



L S |+a 7 , n = 0,±l,±2..., (14) 



j(h,w) 



El + fcn+n_ 2 
2 Sm 

Ufc=±l 



A + + kA. 



1/2 



(15) 



A± = (d/hv)\h±w\, n± = (h±w)/\h±w\ (16) 

which has the standard structure of the equation, for An- 
dreev levelsllj in current-carrying SNS junction!] An ad- 
ditional term (77, where a = ±1 indicates the spin di- 
rection, describes Zeeman splitting of the Andreev levels 
renormalized by the SO interaction. In terms of the BdG 
wave mechanics, the spin effects change phase relations 
between the wave functions of the incident and retro- 
reflected quasiparticles. This produces oscillations of the 
Andreev levels with the change of the interaction param- 
eters h, w which enter the oscillating phase shift 177 in 
Eq. (14). As a result, the Josephson current in Eq. (1) 
also reveals oscillations with h, w as the manifestation of 
complicated interference between partial supercurrents 
carried by individual Andreev bands. In this sense, the 
effect represepts a spin-induced analogue of the Fraun- 
hofer patterntll in a planar Josephson junction in exter- 
nal magnetic field. Note that in the absence of Zeeman 
field (h = 0), all spin effects disappear: 7(0,10) = 0. 

The most striking manifestation of the SO interaction 
itself is the anisotropy of Andreev levels with respect to 
the direction of the Zeeman field h; in the absence of 
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the SO interaction (w = 0), E n ($>) depends only on the 
modulus of h: 

~f(h, 0) = arcsin [sin (hd/Tiv)] . (17) 

It is helpful to consider this anisotropy for a single elec- 
tron state with fixed direction of the SO vector w. As 
follows from Eqs. (14)-(16), the energy of Andreev level 
depends only on the angle between h and w and their 
moduli, being insensitive to the rotation of h around w. 
At h || w, the effects of SO interaction vanish, as in Eq. 
(17). These conclusions can be extended to the angle 
dependence of the Josephson current in Eq. (1) in a nar- 
row 2DEG channel which holds a single electron mode 
(p = pf, P y = 0). In this extreme case, the vectors w 
of all electrons are directed along the y-axis and create a 
fixed reference frame for the Zeeman vector h. 

At arbitrary length of the 2DEG region, Eq. (14) can 
be solved only numerically. Below we consider an analyt- 
ically solvable case of 2DEG channel much shorter than 
the coherence length £o, when the left-hand side of Eq. 
(14) is negligibly small, 

E X {<P) = sAcos($/2 + cr7), (s,a) = ±l. (18) 

At 7 = 0, Eq. (18) describes two spin-degener.ated An- 
dreev levels in a superconducting constrictionta which 
traverse across the whole energy gap with the change 
of $ and intersect each other at $ = 7T.E3 The spin ef- 
fects split each level into two spin-dependent terms and, 
in addition, expand them into four energy bands which 
transfer the Josephson current 

r ._ eA f PF d Pv v-^ [ A /$ V 
m = TU 2^ tanh L 2T C ° S U + ^1 

x sin($/2 + (T7) . (19) 

At h = 0, Eq. (19) is reduced to a 2D analogue of the 
current-phase relationshipEJ for pure constriction. 

The set of curves calculated numerically at T = 
for various directions and magnitudes of Zeeman field 
combined with SO interaction is presented in Fig. 2b-d 
in dimensionless variables 

H = (hv F /d)h, W = {hv F /d)w, (20) 

in comparison with those plotted for W — in Fig. 2a. 
The common features of these dependencies are drastic 
variations of the shape of at H ~ 1, and the rapid 
change of sign of the derivative dl/d& at $ = ir in small 
field H . It is interesting that the curves at W = and 
at W = 1, H || y (Fig. 2a, b) are similar each to other, as 
well as the curves for W = 1, H \\ x and W = 1, H \\ z 
(Fig. 2c, d). This reflects the results of our analysis of the 
anisotropy of Andreev levels in ID case which appears to 
be qualitatively applicable for 2D system: the SO effects 
are relatively small at if || y and approximately isotropic 
under rotation of Zeeman field around the y-axis. 




FIG. 2. Current-phase relationships /($) at T = nor- 
malized on the critical current at H — for various directions 
and magnitudes of the Zeeman field and intensities of the SO 
interaction: H = (1), 0.4 (2), 0.8 (3), 1.2 (4), 1.6 (5). 




FIG. 3. The set of dependencies of the critical current 
I C (H, T)/J c (0, 0) on the dimensionless magnetic field H, at 
W — 1, H \\ x, y, z (solid curves) and W — 0, H directed 
arbitrarily (dashed curves). Upper pairs of curves were cal- 
culated for T = 0, lower pairs - for T = 0.9T C . 

In Fig. 3 we present oscillations of the critical current 
I C (H) depending on the HF field direction and SO inter- 
action intensity, for two values of temperature. At the 
vicinity of T c , where has a single harmonic, 

?reA 2 f PF 

7(*)«/ m sin $, I m = / dp y cos2 7 , (21) 

{Zirnyic j_ PF 

the critical current I C (W, H) = \I„\ turns to zero period- 
ically with H like in SFS system.EI 

The positive sign of dl/d$ at $ = it, which occurs 
at H ^ (Fig. 2), means that this state can be stable 
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and may produce persistent current in the ground state 
of a superconducting loop with high enough inductance 
(7r-junctionEil). On the other hand, the negative sign of 
dl/d& at $ = 0, which occurs within the certain field 
range at W = or at W ^ 0, H || y, signifies instability 
of usual ground state with $ = (note that the SO 
interaction stabilizes this state at H _L y, at least at 
T = 0). The results of a numerical analysis of stability 
of states $ = 0, 7r within the whole temperature range 
are shown in Fig. 4. 



cated from a superconductor with low T c < IK. This re- 
striction can be significantly softened in long (d ^> £o) S- 
2DEG-S junctions where the interaction energies should 
be comparable with the small distance between Andreev 
levels: h,w ~ hvp/d <C A. 



8' 





(a) W=0 


mm 








^^SlllilSI!IIISilil!llilllllilS 



4 


ib! W 1. II y 


3 

l | 










n 









0.5 



1.0 



.0 0.5 
T/T„ 



1.0 



f (c)W=l,H||x 




.0 



0.5 



1.0 



[ (d)W=L Hllz 







0.5 



1.0 



FIG. 4. H-T/T c diagrams of stability of S-2DEG-S phase 
states with $ = and $ = n for various directions and mag- 
nitudes of the Zeeman field and intensities of the SO inter- 
action. Within the gray regions, only usual 0-state is stable: 
d//d$(0) > 0, dl /d4>(7r) < 0; blank regions correspond to the 
stability of the 7r-state only: dl/d&(0) < 0, d//d$(7r) > 0; 
within the dark regions both derivatives are positive. 

In summary, we have shown that the Josephson current 
in a mesoscopic S-2DEG-S structure is highly sensitive to 
the combined action of the Zeeman field and spin-orbit 
interactions. In particular, the critical current reveals 
oscillations and anisotropy with respect to the Zeeman 
field B, and the regions of stability at $ = ir (like in ir- 
junctions) emerge. We assumed hyperfine interaction of 
electrons with polarized nuclei as the source of electron 
spin polarization, though similar effects should be ob- 
served in external magnetic field lying in the 2DEG-plane 
(to avoid orbital effects). In order to access the regime 
of strong interaction (H ~ 1, W ~ 1) in short 2DEG 
bridge (d < £g) considered here, the interaction ener- 
gies h, w of 2DEG should exceed A. Since the HF and 
SO interaction magnitudes in GaAs/AlGaAs heterojunc- 
tions reach at most IK in temperature scale, the banks 
of short S-2DEG-S structure should be preferably fabri- 
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